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Abstract 

The behavior of both the survival S{t) and nonescape P{t) probabihties at long times for the one- 
dimensional free particle system is shown to be closely connected to that of the initial wave packet 
at small momentum. We prove that both S{t) and P{t) asymptotically exhibit the same power-law 
decrease at long times, when the initial wave packet in momentum representation behaves as /c™ 
with m = or 1 at small momentum. On the other hand, if the integer m becomes greater than 
1, S{t) and P{t) decrease in different power-laws at long times. 
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The study of the space-time evolution of the wave packets is very significant for un- 
derstanding of the scattering phenomena and attracting many researchers in the various 
fields. It is then helpful to know the complete information about the free particle system. 
For the one dimensional case, if the Gaussian wave packet is chosen as the initial one, the 
wave packet il){x^t) decreases asymptotically as t'^^"^ at long times. However, it has been 
recently found that the maximum of wave packet does not necessarily behave as t~^/'^ for 
non-Gaussian initial- wave-packet. In fact, a slower decrease than t~^/^ can be found for the 
powe.law tan wave packet Wk, and a faster decease than r^n _ 

n 

packet which vanishes at zero momentum 4]. These facts remind us of a naive question of 
how the characteristics of the initial wave packet affect the long time behavior of the wave 
packet or related quantities such as the survival S{t) and nonescape P{t) probabilities. The 
former stands for the probability of a state still being in the initial state at a later time t. It 
is widely used for the decaying systems (see, e.g., Refs. and references therein). 

The latter is the probability to find a particle in a specific region under consideration at a 
later time t. It is also used to the decaying systems (see, e.g., Refs. Such a 

question for S{t) was answered in the following sense. For a one dimensional free particle 
system, it was shown that S{t) behaves asymptotically like t^^™^"^, when the initial wave 
packet in momentum representation behaves like near the zero momentum = with 

n n 

an arbitrary nonnegative integer m Hence, as is pointed out in Ref. the small- 

momentum behavior of initial wave packet has a crucial role to determine the long-time 
behavior of the survival probability. However, such strict structures for the wave packet and 
the nonescape probability have not been clarified completely. 

In this work, we consider the asymptotic behavior of a free wave packet ip{x,t) at long 



times for the one dimensional case, assuming as in Ref. that the initial wave packet 
behaves like k"^ at small momentum. The asymptotic behavior is evaluated at every position 
X unlike the studies in Refs. |ll Q, 0, Q]. This advantage enables us to discuss whether the 
asymptotic behavior of the wave packet has the position dependence. In addition, we are 
able to calculate explicitly the asymptotic form of not only S{t) but also P{t) at long times 
jl^ . We then examine and clarify the difference between the long time behavior of the S{t) 
and P(t), according to the small- momentum behavior of the initial states. Remark that a 
comparison between the long time behavior of S(t) and P(t) was already made in Ref. [7[ for 
the potential systems in another context, though the analysis therein was correct for S{t), 
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ip{x)ip{x, t)dx 



(1) 



3ut not for P{t). The correct result for P{t) turned out to be the behavior (see, Ref. 

13^ and references therein). 

For the one- dimensional free particle system with the Hamiltonian Hq = 
— {h^/2M)(P/dx'^, we here define the survival probability S{t) of the initial state (wave 
packet) ip as 

S{t) := |(V',e-'*^°/V)r = 

where ip{x,t) = {e~^^^'^/^il>){x) and the bar (") denotes complex conjugate. il){x) is assumed 
to be square integrable. S{t) is the probability that the state at a later time t is found in the 
initial one. We also define the nonescape probability P(t) as the probability that a particle 
initially prepared in the state if) is found in a bounded interval [a, h] on the line at a later 
time 

P{t):= [ mx,t)\'dx. (2) 

J a 

In order to estimate the asymptotic behavior of il!{x,t), we first refer to the explicit 
solution to the Schrodinger equation, 

^{x,t) = -j= / e"'^e-'""''/^^^{k)dk, (3) 

where k = p/h and the ipi^) is the initial state in momentum representation. To see 
the long time behavior of the solution, one can consider the asymptotic expansion of the 
integral in Eq. (jHl), using the phase stationary method [3l as used in Ref. PiQl or making 
an integration by parts for the Fourier integral Then the differential coefficients of V'(fc) 
at k = naturally appear. However, to take account of the x-dependence in the asymptotic 
behavior of ilj{x,t), it may be convenient for us to start with Eq. Indeed, expansion of 
the exponential function in Eq. Q immediately leads to the asymptotic behavior of ilj{x,t) 
with the x-dependence. It reads 



where Gj is the integral operator [l^ defined by 



{Gj^){x) := I \x -yyil){y)dy. (6) 



Here we assume the exchange of the order of summation and integration to be allowed. Note 
that {G2j4'){x) in Eq. (0) can be described in terms of the differential coefficient of ip{k) at 
= as same as the result reached from Eq. Q. This is seen from the following formal 
expansion of i^{k), 

2 /"oo 



OO / . r \ o POO 



(7) 



where ip^'^\0) = i^{0). This implies that 



dki 



k=o v27r J-oo 
Then, we can rewrite the {G2jip){x) in Eq. ® for j = 2j as 



y^^{y)dy. 



(8) 



(G2,^)(X) 



2j 



/27r ^/2j 



2[(2j)!j-V- 



^«^{«)(0)(-x)2j'- 



(9) 



Substituting Eq. © into ((H) and we can obtain the asymptotic behaviors of S{t) and 
P{t) at long times as, 

2 



and 



P{t) 



-i)^--ir(j + i/2), 



(10) 



(11) 



respectively. In Eq. fllUp . {ilj,G2jip) is also described in terms of the differential coefficients 

ip^^\0) as. 



(2?)! ^ \ 

^ -J' n=0 



^^•^U(2,-n)(0)^W(0). 



(12) 



We now consider such a special case that the initial wave packet ip{x) satisfies 



^(A;) = Oik"") as ^ 0, 



(13) 



where m = 1,2,.... We notice from Eq. ((Tj) that the condition (fT^ is equivalent to the 
condition 

^(^■)(0) = 0, for J = 0, 1, . . . , m - 1. (14) 
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Note that the condition ()14p causes the S{t) to behave hke t ^"^ ^. To confirm this assertion, 
it suffices to show that the condition (|T^ imphes the next condition 

(V-, G2j^) = 0, for J = 0, 1, . . . , m - 1, (15) 

and vice versa [ll|. In fact, substitution of Eq. (fT3jl into (fTUj) surely leads to Sit) ~ t"^"^"^. 
We brieffy show the equivalence between the conditions ()14|) and ()15p. The fact that Eq. ()14|) 
implies Eq. fll5|l follows straightforwardly from Eq. ()12j) . Conversely, if Eq. (fTH|l holds, we 
have from {^jj, Goip) = that ip^'^^O) = [see Eq. p2j) ]. Then, we also have from ■iIj^^\0) = 
and ii), G21P) = that ^'^^^(0) = [see Eq. ^ again]. As the same way, we can recursively 
show Eq. (^U), and the proof is completed. Under the condition (fT^ . we see that the first 
non- vanishing term G2mV') the summation in Eq. (fTUI) is reduced to 

{^,G,^ij) = ^ > |^(-)(0)r. (16) 

Then, we obtain the asymptotic behavior for S{t) as |ll| 

Note that this formula is also seen to be valid for m = 0. Let us now examine how the same 
condition (fT^ [or ()14p ] affects the asymptotic behavior of ilj{x,t). Under the condition ()14p. 
we see that Eq. ^ reads 

iG2ji')ix) = for all x G M, (18) 

where j = 0, 1, . . . , m — 1. The m — 1 is the largest integer satisfying 2(m — 1) < m — 1, and 
the m turns out to be 

{m/2 for even m 

(19) 
(m + l)/2 for odd m 

Equation (fTH|) consequently implies that the asymptotic expansion © for the wave packet 
reads 

as t — 00. One also see that this formula is valid for m = 0. By using Eq. (fT^ . the first 
non- vanishing term {G2m4'){x) can be reduced to a simple expression as 



{G2^^)ix) = ^^^^"HO), (21) 



for even m, or 

{G^^)[x) = 2[(^ + i),] [^^"^'H0)+^(^ + W'-HO)], (22) 
for odd m. Substituting Eq. into Q, we can also derive the asymptotic behavior for 

p{ty. 

^^^^ = Jtt/m'^-j ! l(G..^)(x)prf. + 0(t— ), (23) 
as t oo. The above formula is also expressed in terms of the differential coefficients 
?/'(^)(0), by using Eq. (jHl) [or 

It is worth noting that, in the case of m being odd, there is a possibility to find a special 
position, denoted by ^o, where {G2m4'){x) in Eq. ()22|) vanishes. This means that at that 
position the ip{x,t) follows the power-law in the next order. However, this matter may be 
regarded as exceptional, because a point in the entire line has only zero-measure. From Eq. 
(J221), ^0 is given by 

^0 = z^C'n+i)(0)/[(m + 1)V^(")(0)], (24) 

and must be real. We can find such a for the initial wave packet, e.g., 
]Sf^k'^Q-°-o'^(^-^oY I'^-'^^o^ ^ where oq > 0, k^, xq G M, and Nm being the normalization constant. 
In this case, the .^o is given by XQ + iaQ^kQ. Then, it becomes real if and only if k^ = 0, which 
leads to ^0 = 2^0, the center of the initial wave packet. Note that such a special position in 
Eq. ()24|) if any does not have an influence on the asymptotic form of P{t), because P{t) is 
obtained by the integral of \{G2m'ip){x)\'^ ■ 

Let us now consider and compare the long time behavior of S{t) and P{t). We see from 
Eq. (|19p that m and m are different when m > 2, and this fact directly affects the long 
time behavior of S{t) and P{t). When the initial state ip satisfies ipi^) = 0{k"^) with an 
arbitrary integer m > 2, S{t) goes asymptotically like t~^™~^, whereas P{t) does like 
for even m or like for odd m. For a large m, S{t) decreases much faster than P{t). 

We also see that m = m + 1 for odd m. This means that, in the case of an odd integer m, 
unlike S{t), P{t) decreases in the same power law under both the conditions, i^i^k) = 0{k"^) 
and ^{k) = 0(fc™+^). 

To illustrate the difference in the long time behavior of S{t) and P{t), we choose the 
three initial wave functions (po{x), 4>i{x), and (p2{x), defined by 

^^(k) = iV^fc™e-'^«'('=-'=«)'/2-«ofc^ m = 0, 1, 2. (25) 
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FIG. 1: (a) S{T) and P{T) (square and diamond, respectively) of the wave function (p^ in Eq. ()25() . 
and their asymptotes predicted by Eqs. (|17j) and (sohd and dashed hues, respectively), where 
T = ht/2Mao^ being the reduced time. In this case, S{T) and P{T) show the same power decay 
behavior like at long times, (b) S{T) and P{T) of the wave function 01, and their asymptotes. 
The notations and symbols are the same as those in (a). S{T) and P{T) exhibit the same power 
decay, however they behave like T^^ instead of T~^. Here we set /cq = 0.0 and xq = 0.0, in (pQ and 
01, and a/ao = —2.0 and b/ao = 2.0 in P(T)'s . 

These are the same ones considered after Eq. ()24j) . They behave like 4>m{k) = 0{k"^) for 
small k. Figure n shows the time evolution of S{t) and P{t), and their asymptotic forms 
predicted by Eqs. (fTTj) and (^Hj). The initial states 0o and 0i are used in Fig. [T] (a) and 
(b), respectively. It is clearly seen that in Fig. [T] (a) S(t) and P(t) behave asymptotically 
like at long times, and in Fig. ^ (b) like t^^. In these cases, the difference between 
the behavior of S{t) and P{t) is not found. On the other hand, we notice that in Fig. 
El S{t) and P{t) for the initial state 02 differ asymptotically at long times. The former 
behaves asymptotically like t~^, however the latter behaves like t~^. In our calculation, we 
have chosen a set of parameters ko = 0.0 and Xq = 0.0 for the three initial states, and 
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FIG. 2: S{T) and P{T) of the wave function (/>2, and their asymptotes, where T = ht/lMao^ . The 
same notations and symbols as in Fig. ^ are used. In this case, S{T) and P{T) exhibit different 
power decays at long times. The former behaves like while the latter like T~^. Here we set 
/co = 0.0, and xq = 0.0, in (j)2, and a/ao = —2.0 and 6/ao = 2.0 in P(T). 

a/do = —2.0 and b/ao = 2.0 for the interval [a, b] for P(t). Then, as is seen from Figs, ^and 
El -P(O) ~ 1, i.e., the initial states are well localized in the interval. 

In conclusion, we have considered for every position the long time behavior of the wave 
packet moving freely in one dimension, according to the characteristics of the initial wave 
packet at small momentum. We then have found that the asymptotic power of t obeyed by 
the wave packet is constant everywhere, at most excluding one position ^o- We also have 
obtained the asymptotic behavior of the nonescape probability at long times, and compared 
that of the survival and nonescape probabilities. It is of interest that they can decrease in 
the different power laws depending on the initial states, in spite of the apparent similarity 
between their physical meanings. Our derivation can be easily extended to an arbitrary 
dimension, by starting with Eq. (0} in a corresponding dimension. In these analyses, we 
assume that the exchange of the order of summation and integration is admitted in the 
formal expansions in Eqs. (0) and (|7j). Indeed, this assumption can be rigorously guaranteed, 
when we make the same discussion with the finite series involving an appropriate remainder, 
instead of Eqs. © and ((Zj). In any such procedure to keep the validity of the formula, e.g., 
()20|) . what should be satisfied at least is that all of the differential coefficient ^^^\0) with 
j up to m (or m + 1) is finite for even m (or odd m). See Eq. (PT|) [or ((221)] • It should 
be noted that this condition also implicitly implies that limk^+o %()^^\k) = limk^^Q i/j^^^k) 
for j = 0, 1, . . . , m (or m + 1). These conditions are satisfied by those ip^s who are rapidly 
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decreasing functions as in Eq. ()25|). However, such a circumstance is not always valid for an 
arbitrary initial wave packet, e.g., the wave packet with the power-law tail [ll O, or that 

causes |V^(0)| = oo and the latter does linifc^+o ^^'■'"-'(/i;) 7^ 
limfe_^_o 'i/'^™'* (A;). It is then significant to consider how our results are modified for such initial 
wave packets. Furthermore, it is important to extend our consideration to the potential 
systems. In particular, it is relevant to examine in that case the possible influence of the 
characteristics of the initial states in the long time behavior of the survival and nonescape 
probabilities. In fact, such an attempt has not been done in previous investigations. An 
extension may be realized by starting, instead of Eq. (jSl), with the asymptotic expansion 
of the wave packet at long times for the short-range potential systems, attained by several 
methods (see, for example, Refs. 



17| and references therein). 
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